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RESUMEN 

Se aplico el metodo de vectores principales de disco y extincion a los datos 
fotometricos en el cercano infrarrojo de los cumulos estelares jovenes de las 
regiones Tauro-Auriga y Orion. Bajo la suposicion que la edad del cumulo 
esta representada por la mediana de la distribution de edades se obtiene una 
estimation de las masas estelares individuales. La funcion initial de masa 
(FIM) obtenida para estos cumulos exhibe una gran similitud con la obtenida 
por metodos espectroscopicos y fotometricos y puede constituir una repre- 
sentation robusta de la FIM. El metodo permite encontrar la contribution 
por extincion y disco para cada estrella. La distribution global de la ex- 
tincion para el cumulo de Orion se compara bien con resultados previos. Se 
encuentra un dominio en la frecuencia de estrellas T Tauri con disco. 

ABSTRACT 

We applied the extinction-disk-principal vectors approach to near infrared 
photometric data of the Taurus- Auriga region and Orion Nebula young stel- 
lar clusters. By assuming that the cluster age is represented by the median 
value of the age distribution we are able to derive the distribution of stellar 
masses. We showed that the resulting initial mass function (IMF) for these 
two young stellar clusters compares remarkably well and might be a robust 
representation of the IMF obtained by spectroscopic or photometric meth- 
ods. The method also yields extinction and disk contribution for each star. 
The overall extinction distribution for the Orion cluster is analyzed and com- 
pares well with previous work. The frequency of T Tauri stars with disks is 
dominant. 

Key Words: Stars: formation — Stars: fundamental parameters — Stars: 
low-mass — Stars: pre-main sequence 

1. INTRODUCTION 

Masses of pre-main sequence (PMS) objects are quite difficult to deter- 
mine, and the only direct way to assess them is throug h the analysis of their 



dynamical parameters. Following I Cohen fc Kuhil (jl979l ) locating objects in a 



temperature-luminosity H-R diagram along with evolutionary models yields 
masses. Since stellar temperature is estimated from colors, and these are 
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affected by circum stellar disks, PMS te mperatures are best estimated spec- 
troscopically (e.g., Luhman et al. 2005h . 



But PMS objects are usually faint 
and found in dusty environments and obtaining their spectra becomes quite 
difficult. This has lead many authors to explore the initial mass function 
(IMF) of young stell ar clusters in the in frared, via the luminosity function 
in the i-T-band (e.g.. iMuench et al.ll2002l ). assuming that it reproduces the 



true IMF of a cluster. Color-color (CC) and color-magnitude (CM) near- 
infrared diagrams of young stellar clusters show that intrinsic near-infrared 



excess emission (Lada & Adams 1992 


1997; Hillenbrand & Carpenter 


2000). 


Lopez-Chico & Salas ( 


2007), hereoi 



stars can be obtained using JHK photometry and both CC and CM diagrams. 
This is the result of the analysis of two principal vectors, one produced by the 
disk excess, D and the other by interstellar extinction, X, if stellar ages are 
known and a particular set of evolutionary tracks is assumed. 

In this paper we continue the development of a new approach in determin- 
ing masses of pre-main sequence stars from near-infrared photometry. Our 
goal is to strengthen the results reported previously in LS07. First, we refine 
the values of disk excess coefficients given in LS07 by showing that the method 
can be extended to the I and L filters, and that these coefficients scale well 
with wavelength. This analysis is presented in section 2. Second, we test 
whether the method can be used to extract the Initial Mass Function (IMF) 
of young stellar clusters from infrared photometric data alone. Sect. 3 shows 
that if the age of a young cluster is known the distribution of stellar masses 
can be obtained from their JHK photometry and a set of PMS evolutionary 
tracks. As a proof of this statement we apply the method to the well studied 
Taurus-Auriga region and the Orion Nebula Cluster. It is shown that the 
median age of the clusters produces an excellent agreement with previously 
known IMFs. These results are reported in Sects. 4 and 5, while a summary 
of our results is presented in Sect. 6. 



2. METHOD DESCRIPTION 



Lopez-Chico fc SalasI ( 2007 ) showed that masses of T Tauri stars can be 



obtained using their JHK photometry and their location in both CC and 
CM diagrams. This is the result of the analysis of two principal vectors, one 
produced by the disk contribution, D, and the other by interstellar extinction, 
X. LS07 show that the vector X can be defined as the extinction vector 



corresponding to ^4y = 10 with components given in lRieke &: Lebofskvl (|1985l ). 
While the vecto r D is obtained via J, H, and K magnitudes obtained from 
D'Alessio et al.1 (|2005f ) models of accretion disks irradiated by a central T 
Tauri star. These base vectors were thus given the following components in 
CM ( K vs. (J - K)) and CC ((J - H) vs. (H - K)) diagrams: 



D cm = (1.014, -1.105) ; X cm = (1.68, 1.16) 
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TABLE 1 

COEFFICIENTS X x AND Y x FOR 
EXTINCTION AND DISK 
CONTRIBUTIONS 



Filter 


x£ 


yx 


Vxl 


I 


4.82 


0.064 


0.08 


J 


2.82 


0.15^ 


0.15 


H 


1.75 


0.40 c 


0.43 


K 


1.12 


1.19 c 


1.12 


L 


0.58 


4.86 


4.93 



"From Rickc & Lebofakvl jl98Sft 

''Using y\ tx A 3 . 

c From lL6pez-Chico fc Salad {200' 



(1) 

DZ = (0.767, 0.221) ; JQ C = (0.61, 1.07) 

The position of a particular star in these CC and CM diagrams becomes then 
a vector sum in each diagram that can be reduced to a closure relation plus 
the following set of linear equations: 

J = J + axj - (3y } 
H = H + as H - /3y H (2) 
K = Ko + ax K - f3y K 

where a and (5 represent the amount of extinction and infrared excess con- 
tributions to the individual magnitudes, x\ and y x are the corresponding 
extinction and disk coefficients for each wavelength A, and the sub-indices 
denote the intrinsic photospheric magnitudes. Values for x\ and yx are ob- 
tained directly from the components of base vectors given in eq. (1) and are 
presented in Table 1. 

An advantage of presenting the equations as a set of linear equations is 
that it can be easily explored with any set of three photometric filters at a 
time, that is, we may use filters UK or JKL instead of JHK, provided the 
corresponding closure condition is fulfilled and that the excess colors behave 
as vectors. We found that this procedure is also feasible for filters I and L in 
addition to the JHK set, and so the numerical values of the yx coefficients 
for I and L are included in Table 1 . We show in Figure Q] that the relation 
between log yx and log A can be well represented by a power-law index equal 
to 3. This close relation gives us confidence to refine the yx values as given in 
the last column of Table 1, which are the values that we used in our analysis. 
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Fig. 1. Plot of log yx as a function of log A shows that for IJHKL the data are well 
represented by a power-law of index=3, shown as a solid line. 



As is shown in LS07 stellar masses can be obtained from the solutions of 
the minimization of the quadratic error Err given by 

Err 2 = ^ I (3) 

A=r 

where m° bs is the observed magnitude at each wavelength A and m\ is the 
magnitude that should be observed according to the model: 

m\ = m° (mass, age) + d + ax\ — (3y\. (4) 

In eq. (4) the stellar magnitude is the sum of the absolute magnitude for a 
certain mass and age, plus the distance modulus d to the object, plus the ex- 
tinction correction and the infrared excess contribution from the circumstellar 
disk. The error is then minimized with respect to its four parameters: mass, 
age, a and (3. This minimization procedure is a set of linear equations for a 
and /3. To deal with mass and age we compute Err for each single mass and 
age taken from pre- main sequence evolutionary models. Then one seeks the 
minimum of Err consistent with the additional constraints a > and f3 > 
obtained for each star. This method is called the extinction-disk principal 
vector, XDPV method. 



To test the XDPV method we will use only the PMS tracks from lD'Antona fc Mazzitelli 



(|1997l ) and http : //www .mporzio . astro . it/~dantona/prems .html that pro- 



vide luminosities and temperatures for pre-main sequence stars in a wide low- 
mass range, from 0.017 to 3 M and ages from 10 4 to 10 s yr. The use of other 
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evolutionary tracks, (e.g. Palla fc Stahler 19991 ). has been discussed in LS07. 
These authors argue that these tracks are found to produce similar mass re- 
sults within 20%. Luminosity and effective temperature of the evolutionary 
tracks arc converted into absolute magnitudes using b olometric corrections 
and normal colors given in lKenvon fc Hartmannl (jl995l ). 



3. MASSES OF YOUNG STELLAR CLUSTER STARS 

In LS07 we showed that if the masses of young stars are known, e.g. spec- 
troscopically, the proposed XDPV method is a powerful tool to derive their 
ages or conversely, if the ages are known the masses can be extracted. How- 
ever, both cannot be derived simultaneously. This is due, as pointed out in 
LS07, to the fact that the minimal of Err as a function of mass and age is a 
region that resembles a long and narrow canyon (c.f. their Fig. 7) that spans 
a wide range in masses and ages, providing a continuum of possible solutions. 
It is then necessary to specify an age for the PMS object to determine its 
mass. 

However, when dealing with a cluster of stars with similar ages, a repre- 
sentative age of the cluster may be used for all the individual members. Doing 
so requires a compromise age that would compensate the errors in mass by 
assuming a too young age for some objects with those errors derived from as- 
suming an older age. The median age bisects the age histogram in equal parts, 
so that an equivalent number of members is either younger or older. Further- 
more, the median is a robust indicator of the central tendency, its value is the 
same whether one uses the histogram of ages in linear or logarithmic values, 
and in general is a better choice when we only have one number to specify a 
distribution. For these reasons, we have chosen to use the median age of each 
cluster. 

We will show that this age selection gives consistent results for the resulting 
mass histograms of the well studied young stellar clusters of Taurus- Auriga 
and Orion, through the comparison of the IMF obtained by other authors 
using spectroscopic or photometric methods, and the IMF obtained applying 
our algorithm. We note that in the young stellar clusters studied here, the 
median age has a lower value than the mean, due to a fraction of younger stars 
in the tail of the distribution. For Taurus- Auriga we obtain log median = 5.8 
(0.63 Myr) and log mean = 5.95 (0.89 Myr), while for Orion log median = 5.6 
(0.4 Myr) and logmean =5.8 (0.63 Myr), as is described in detail below. 
These values are significantly younger than what is commonly assumed as the 
age of the clusters. But this has to be that way, because one usually refers to 
the epoch at which star formation began, which is the upper limit of the age 
distribution. 

The method produces very good results as is shown below. Although it 
may seem a disadvantage to have to specify a cluster age, it is still far less 
information that the requirement of spectral types in spectroscopic methods, 
or a parametric description of the age distribution required in some other 
photometric methods. 
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4. TAURUS-AURIGA 



As pointed out by lKenvon fc Hartmannl () 19951 ) . hereon KH95, the Taurus- 
Auriga is an ideal laboratory to study low-mass star formation, dominated by 
low-mass stars with little extinction, and so has been the subject of many in- 
vestig ations. The list of known members of the region has grown through the 
years JCohen fc Kuhilll979[ iHerbig fc Bel3ll988l Il995l iKenvon fc Hartmann 



1995 



Briceno et al 



2006; Scclsi et al 



2006 



Guieu et al 



20021 ; iLuhmanl 12004 iLuhman et al. 

with an increasing emphasis towards completeness. 



As we discussed above, to estimate masses the XDPV method requires 
age estimates in addition to near-infrared photometry. We will take this age 
from the histogram of ages presented by KH95 (their Fig. 16). They have 
derived these ages together with masses from a full set of visible and in- 
frared photometry and spectral types for known Taurus- Auriga members in 
1995. Their age distribution spans from log(age(yr)) = logT = 4 to 6.5 
and we calculate a mean value of 5.95 and a median value of 5.8 (T=6.3xl0 5 
yr), which on one hand is smaller than the commonly assumed value of 6.3, 
and also smaller than the mean. We may go on to compare masses derived 
by our method. KH97 present spectral types and effective temperatures for 
139 stars. They use 103 of them to construc t a mass histogram from H- 
R diagrams with D'Antona fc Mazzitellil |l994h tracks with CMA opacities. 
We took 189 stars with JHK photometry from their list, and calculated 
masses as described in Sect. 2 above, using the median age derived (logT= 5.8) 
from KH95, and the evolutionary models of lD'Antona fc Mazzitelli (1997) and 
http : //www . mporzio . astro . it/~dantona/prems . html 

We compare mass histograms in Fig. [2] To do so, we have re-binned 
KH95's histogram in logarithmic 0.3 dex bins, by assigning random masses 
to each star within its own bin, and then re-binning in the logarithmic bins. 
We repeated this process 100 times to produce the mean histogram that is 
shown by dotted lines in Fig. [51 while ou r mass histogram is shown as solid 
lines. For reference the Miller-Scalo IMF ( Miller fc Scalolll979h is shown as a 
dashed line. The general shape from 0.3 to 1 Mq is remarkably similar in both 
histograms and both peak at the same value where a turnover is observed. The 
excess number of stars in our analysis (189 of 103) appear in three regions 
of the distribution. Some are in equal proportions in the three bins from 
logAf =0.1 to -0.5 Mq, some lie in the high-mass end of the distr i butio n 
which nevertheless agrees quite well with the IMF of iMiller fc Scale! |l979h . 
and the majority are distributed in low-mass bins log M /Mq < -0.8, that most 
probably were too faint to provide a reliable spectral type by KH95. However, 
this low-mass region is not unbiased or complete, as is pointed out in their 
paper. 



Briceno et al.1 (|2002l ) ; ILuhmanl (|2004h and ILuhman et al.l (|2006j ) have paid 



attention to this fact, and have conducted spectroscopic surveys of several re- 
gions in the Taurus- Auriga clouds in order to identify all low-mass stars (many 
of which are brown dwarfs) to get complete unbiased samples of the IMF. 
This IMF is the result of photometry and spectroscopy to determine lumi- 
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Fig. 2. IMF in Taurus- Auriga from lKenvon fc Hartmann ( 1995) data processed by 
the XDPV method (shown as solid lines) compared to lKenvon fc Hartmann! 1 1995h 
mass histograms (dotted lines). 
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Fig. 3. IMF in Taurus- Auriga from iLuhman et al.l (12006 ) data as obtained by the 
XDPV method (shown as solid lines) compared to iLuhmanl 1 2004 ) 's IMF (dotted 
lines) . 



nosities and effective temperatures, complemented with Baraffe et al. ( 1998h 
evolutionary tra cks to der i ve ma sses. The more consolidated example of this 
I MF is given in ILuhmanl |20o4 . We have taken 2-MASS JHK data given 



Luhman et al 



1 2006) for one of such complete regions that are consis- 



tent with ILuhmanl (|2004j ) IMF. In this new article, Luhman incorporated 
some 20 newly discovered brown dwarfs into a large region that encompasses 
about half of the known Taurus population, and thus constitutes a signif- 
icant sample. With our method we were able to obtain solutions for 125 
objects out of the presented list of 156, again using the median age d erived 
from KH95, and the evolutionary models of iD'Antona &: Mazzitellil ( 1997 ) 
and http : / / www . mporzio . astro . it/~dantona/ prems . html. 

As can be obse rved in Fig. [3] the resemblance between the IMF obtained 
by ILuhmanl (|2004f ) and the IMF obtained using the XDPV method for the 
Taurus- Auriga region is remarkable, with a probability of over 60% of being 
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random samples of the same population. For a Kolmogorov-Smirnov test of 
this result and others, see section 6. 

Furthermore, the XDPV solutions for both data sets ( Kenvon fc Hartmannl 
19951 lLuhmar]|2004l ) presented in Figures 2 and 3, resemble one another better 
that those presented in the original papers. We conclude that this agreement 
strengthens our results. 



5. ORION NEBULA CLUSTER 



5.1. About the IMF 

The Orion Nebula Cluster (ONC) centered on the Trapezium OB stars is 
the richest of any nearby clusters and has been studied extensively. Numerous 
studies have targeted th is cluster to determine its underlying population and 



the a ssociated IMF (e.g. lHillenbrandll997uLuhman et al 



•l ying p 
Jl2000j; 



Muench et al 



20021 ). We will compare our own results to those of two important studies, 
on e spectroscop i c and one based on near-infrared photometry. 



Hillenbrand! (|1997l ) has obtained spectral types for 934 visible stars. This 



information, supplemented by optical photometry, allowed her to populate 
an H-R diagram with pre-main sequence evolutionary tracks, and to ex- 
tract mass and age information. She notes, however, that large uncertain- 

ties arise from the choice o f a particular set of evolut i onary tracks (see also 

Hillenbrand fc WhiteM2004h . Nevertheless JHillenbrandl |l997t ) chooses lLVAntona fc Mazzitelli 
( 1994h evolutionary models to display the ONC's IMF (her Fig. 17), and it 
has become a seminal reference for this region. The derived distribution of 
stellar ages of the ONC population conform a distribution that spans a wide 
range of ages. It starts at logT =3.5 (3000 yr) and increases gradually until 
logT =6.3 (2 Myr), then decreases abruptly and continues at a low con- 
stant pace up to logT = 7.8 (63 Myr). This distribution has a mean value of 
log < T > = 5.84 (0.7 Myr) although some authors quote 0.8 Myr. It has also 
been represented by a constant rate from logT =5 (0.1 Myr) to logT =6 (1 
Myr). In our treatment of the OMC cluster we will choose the median age, 
logT = 5.6 (0.4 Myr), which we believe is the best com promise as is di s cussed 
above for Taurus- Auriga. This is in agreement with Luhman et al. ( 2000l ) 
wo rk that also quotes a median age of 0.4 Myr. 

Muench et all (fcooch developed a Monte Carlo method to model the IMF 



based in obtaining the KLF from a series of probability distributions: extinc- 
tion, infrared excess, age and the IMF modeled as a seri es of power laws. The y 
applied this method to JHK observations of the ONC in lMuench et al.l (|2002t ) . 
The age distribution was chosen as a uniform distribution from 0.2xl0 6 to 
1.4xl0 6 yr. The extinction distribution was derived from (J — H) vs. (H — K) 
diagram by de-red dening sources down to the classical T Tauri star locus of 
Mever et al.l (jl997l ). After this, the infrared excess distribution was obtained 
from the remaining excess in the (H — K) color after subtracting a histogram 
of (H — K) colors of field stars, and this excess color was assumed to arise 
exclusively from excess disk emission at K . From an original set of ^1000 
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sources, they took an extinction limited sample (Ay < 17) of 583 stars, which 

is said to be complete down to 0.017 Mq. 

We took JHK photometry from the iMuench et all J2002) published list, 
assu med a median age of log ? 1 = 5-6 from iHillenbrandl (|l997l ) age histogram, 



usedD'Antona & Mazzitell il (|!997h and http : / / www . mp orzio . astro it/~d antona/ prems . html 

evolutionary tracks, and assumed a distance of 400 pc ( Muench et al. 20021 ). 
With these ingredients we applied our XDPV method to the 699 object where 
no confusion flags are found and photometry is available in all JHK filters. 
We were able to obtain solutions consistent with both a > and (3 > and 
an acceptable Err < 0.3 mag for 612 of these stars. Objects for which no 
solution was found are, for ex ample, stars more massive than those present in 
D'Antona fe Mazzitellil ( 19971 ) and|http : //www .mporzio . astro . it/~dan.tona/prems .html 



evolutionary tracks, maximum of 3 Mq, which includes all the Trapezium 
stars, BN and 6 2 Ori A. From the 612 PMS objects we then selected 578 with 
Ay < 17 to display i n Fig . [4] This compares our IMF ( solid histogram) with 
that of Hillenbrand (1997), points with err or bars, and Muench et al. ( 20021 ) 
shown as a dashed line. The agreement with Muench et al. ( 20021 ) is"excellent. 
that is, the initial slope for high-mass stars, the flattening and the position 
of the turnover, followed by the negative slope in the subsolar mass spectrum 
all agree quite well. The only exception is for the low-mass secondary peak in 
the substellar region (log M/Mq — -1.8). Unfortunately, thos e objects come 



from t he low-brightness peaks in the H and K histograms in IMuench et al 



(|2002f ) but are too faint and red, and therefore, absent in J. Consequently, 
they do not appear in our d ata. 

The agreement with the iHillenbrand ( 1997 ) IMF is quite good in the -0.9 
< log M/Mq < 0.45 range. Massive stars (M>3 M Q ) are missing in our his- 
togram as mentioned above, due to the Unfi tted mass range in t he evolutionary 
tracks. In the log(M/M Q ) < -1 range, the Hillenbrand! ( 1997) survey is mos t 
likely incomplete for sources with Ay > 2.5 |ffillenbrand fe Carpentei]|2000l) . 
The position of the turnover also agrees, althoug h the exact position may 
be at slightly lower masses as has been revised in IHillenbrand fc Carpenter 
(|2000j ) by using updated evolutionary tracks and transformations. 

We conclude that the IMF we obtained with the XDPV method is a ro- 
bust representation of that obtained by other methods. For a more through 
comparison of these distributions see section 6. 



5.2. Extinction and infrared excess 

In addition to the mass data, our method also gives information about 
the extinction a and infrared excess f3 of the sources. We display these two 
quantities as hi stograms in Fi g. and c ompare them with extinctions and 
excesses used in IMuench et al. ( 2002 ) and Hillenbrand! (1997). 

The bottom left panel shows the histogram of extinction a as a solid 
histogram. This quantity can be directly understood in terms of the visual 
extinction, a — Ay /10. It is then easy to compa r e it w ith the extinction 
probability distribution presented in IMuench et alj (|2002l ). shown here as a 
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Fig. 4. ONC's IMF (solid histogram ) compared to 
error bars) and iMuench et al.l (2002) (dashed line). 



Hillenbrand 



( 19971 ) (points with 
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solid line, which is very close to our result. In order to c ompare with the spec- 
troscopically derived extinctions in iHillenbrandl ( 1997 ). we show the dashed 
histog ram for thos e sour ces that were analyzed by us and that are also part 
of the IHillenbrandl (|l997| ) survey, and compare it to the dashed line obtained 
from extinctions derived by her. A di screpancy is notori ous for the very low 
extinction (Ay < 2) sources, where Hillenbrand ( 19971 ) finds most sources 
and our histogram turns over. 

In the bottom right panel we show the h i stogra m for (3 compared to the 



infrared excess distribution in iMuench et alJ (|2002 ). However, there are two 



possibilities for this comparison. If the abscissa in IMuench et all |2002f) Fig. 8b 
is the excess in H — K, then from eq. (2) [3 = -Eh-kAj/k — Va), while if the 
abscissa is taken as the excess in K alone, then f3 = —Ek/uk, and given the 
values in Table 1, they are not equal. We show both possibilities in Fig. 5d. 
The infrared excess is best represented by (3 and we find a better agreement 
when the infrared excess distribution is represented by —Ek, rather than 
Eh-k- 

In the top panels of Fig. [5] we show a and (3 as functions of mass log M. As 
a general rule, no dependence of these two quantities with mass is found, as a 
and (3 acquire all their values for any mass in the range from 0.017 to 3.0 M Q . 
The one exception is observed in the case of a (Fig. 5a): there are no low- 
mass stars at hig h extinctions . This is more notorious in the case of those stars 
also observed by IHillenbrandl (|1997 ). marked as solid dots, where a diagonal 
line in the upper rightmost part of the diagram can easily be drawn. This 
is expected to be the case, since low-mass stars cannot be observed at high 
extinctions for given observation times, and since this fact was not introduced 
a priori in the method, it constitutes another confirming result of the XDPV 
approach. In this diagram it is then possible to extract an extinction limited 
sample, as a rectangular box that l ies below the diago nal line corresponding 
to each sample. It can be seen that IHillenbrandl (|1997l ) survey is unbiased for 
logM/Mfp >-!.! and Ay < 3. in close agreement with her findings (see also 
Hillenbrand fc Carpenter! I2000T ) . On the other hand , the completeness limit 



(200 



2) survey would go up to 



in all three JHK filters in the IMuench et al 
log M/M >-1.8 and Ay <6. Finally, given that the infrared excess (3 is 
due to the disk contribution we find that 300 stars out of 612 (close to 50%) 
presumably possess associated disks {(3 > 0), independently of their mass. 
This is consistent with recent Spitzer studies that show t hat the percentage 



of low-mass PMS stars with disk in ONC is about 50% ( Rebull et al 



200 



6; 



Cieza fc Baliberll20071 ) 



6. STATISTICAL SIGNIFICANCE OF IMFS 

For two young stellar clusters a general concern is whether the individual 
IMFs are different or not, given that the shapes of the IMFs look similar. 
To quantify this we performed a Kolmogorov-Smirnov (K-S) test. First, we 
compared the mostly accepted mass distributions for Taurus and Orion, to 
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log M/M Q log M/M, 




Fig. 5. Top panels: a and as function of mass. In panel a) sources also in the 
survey of iHillenbrandl (|1997l ) are marked as dots. Bottom panels: histograms of 
a and 13 in solid lines, are compared to c) the extinction pro bability distr i bution 
of iMuench et all |2002h solid line, and the extinction data in IHillenbrandl |1997) 
in dashed line to its corre sponding dashed line histogram, d) The infrared excess 
probability distribution of IMuench et al. if considered as E(H — K) (dashed 

line) or if- excess (solid line). 
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show that they are significantly different. Second, we compared the IMFs 
obtained by the XDPV method to both the Taurus and Orion cases. 

To perform the K-S test we developed a routine that generates a random 
population of stellar masses from any given mass distribution iV(logM) using 
the well known accept-reject algorithm. For each distribution that we tested, 
we generated a number of random stellar masses equal to the number of 
stars that were used to obtain the original distribution, in order to keep the 
same stochastic variability that would be expected in a different realization 
of the distribution. For each comparison between different distributions, we 
took 100 realizations of each distribution being compared, and computed 
the two distribution K-S probabilities that the samples are drawn from the 
same parent distributions (null hypothesis). We report the mean of these 100 
comparisons (K-S MP). 

We first tested the IMFs reported by 



(hereon TL) to those of lHillenbrandl (|l997l ) and lMuench et al.l (|2002h . hereon 



Luhmanl (120041) for Tauru s-Auriga 



OH and OM respectively, for Orion. The OH distribution has the advantage 
of being derived spectroscopically (as well as the TL one) , but is incomplete 
for masses lower than log(M/M Q ) < —0.9 as we already mentioned. For 
this reason we limited this distribution to the medium mass range in the 
comparisons. We therefore label it as OHM. Similarly as we already discussed, 
the OM distribution has a secondary peak at the very low mass end which we 
cannot test because the sources involved are below the J filter detection limit. 
Therefore the OM distribution is limited to log(Af/Af Q ) > — 1.6 for all tests 
here. The comparisons between TL and OHM yield a K-S mean probability 
(K-S MP) of 3xl0~ n of being drawn from the same parent distribution, 
while TL and OM gives K-S MP of lxl0~ 4 . These very small values makes 
us confiden t that our test is consistent with previously known results (e.g. 



Lada et al.l (|2008l )) that the distributions are indeed different. 



We then compared the distribution of XDPV masses obtained for the 
Taurus- Auriga region (hereon TX) to TL, obtaining a K-S MP of 21% there- 
fore confirming the posibility of both distributions being equivalent represen- 
tations. The comparison of TX to OHM rules out the null hypothesis with 
a K-S MP of 6xl0~ 5 . It should be mentioned however, that TX is not that 
different from OM (K-S MP of 1%) mainly because in the low mass range 
-1.6 < log(M/M Q ) < -0.9 TL is quite similar to OM (K-S MP = 67%). 
That is to say, the low mass slope of Taurus and Orion IMFs are similar. 
Nevertheless, the agreement between TX and TL is sustained. 

Finally, the distribution of XDPV masses obtained for the Orion region 
(hereon OX) compares with OHM with a K-S MP of 3% in the medium mass 
range, and a somewhat better agreement is seen to OM with K-S MP = 9% in 
the whole range of masses log(M/M©) > —1.6 . On the other hand, OX and 
TL cannot be accepted as a match with a very low K-S MP of 2xl0 -7 . This 
proofs that OX is not consitent with the Taurus IMF but has a fair probability 
of being drawn from the same parent population of the Orion Nebula Cluster. 
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7. SUMMARY 



We used the extinction-disk-principa l vectors approach, wh ich is called the 
XDPV method, reported previously in Lopez-Chico &: Salasl ( 2007 ) to show 
that it is a powerful tool to estimate masses of pre-main sequence stars in 
clusters. The method requires a minimum of information, using as little as 
JHK (or UK or JKL) near-infr ared photometry, suppleme nted by a set 
of PMS evolutionary tracks (e.g. iD'Antona fc Mazzitellilll997l ) and the me- 
dian age of the cluster. For each star in the cluster, we are able to estimate 
the contribution of the extinction vectors A cc in the color-color diagram and 
A cm in the color-magnitude diagram, and the disk vectors D cc in the color- 
color d iagram and D cm in the color-magnitude diagram, using D'Alessio et al 



(2005() accretion disk models grid of spectral energy distributions. The ob- 
served absolute magnitude at each wavelength A of a PMS object is obtained 
via m\ — mP x (mass, age) + d + ax\ — /3y\, where the first term corre- 
sponds to the absolute magnitude of the naked object for a certain mass 
and age, followed by the distance modulus d, the extinction correction (a^) 
and the infrared excess contribution {—(3y\) from the circumstellar disk. It 
is shown that if a representative age of the cluster, such as the median, is 
known, the masses of each individual star can be statistically obtained from 
its near-infrared photometry alone. The XDPV method is tested in the well 
studied regions of Taurus- Auriga and the Orion Nebula Cluster by extract- 
ing their Initial Mass Function. These IMF are in excell ent agreement ( K-S 
test) to tho se given by Kenvon fc Hartmann ( 1995h and Luhmanl (2004) for 



Taurus and Hillenbran 



"dl |l997h 



and 



Muciicl) <•! ill. i 2002)" for Orion. Since 



our algorithm also yields the extinction and disk contribution for each star 
the distributions can be obtained. The overall extinction distribution for the 
Orion c luster is analyzed and compares well with previous work, the compar- 



ison to 
vector is 



Muench et al.l (|2002h shows that the parameter a for the extinction 



Ek rather than k- The frequency of PMS low-mass stars with 
disks, represented by the parameter (3, is about 50% in the Orion sample. It 
is also seen that the number of sources observed decreases for high values of 
extinction a, confirming a well known and expected observational effect, and 
allowing us to draw a complete extinction limited sample a posteriori. We 
conclude that our XDPV algorithm can be applied to study the IMF of young 
stellar clusters, as well as the distributions of extinction and disk infrared 
excess. 



We are grateful to an anonymous referee for very valuable comments that 
aided in the final version of the manuscript. 
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